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Abstract 

A spectroscopic application of the atom laser is suggested. The spectroscopy termed 2PACC 
employs the coherent properties of matter-waves from a two pulse atom laser. These waves are 
employed to control a gas-surface chemical recombination reaction. The method is demonstrated 
for an Eley-Rideal reaction of a hydrogen or alkali atom-laser pulse where the surface target is an 
adsorbed hydrogen atom. The reaction yields either a hydrogen or alkali hydride molecule. The 
desorbed gas phase molecular yield and its internal state is shown to be controlled by the time and 
phase delay between two atom-laser pulses. The calculation is based on solving the time-dependent 
Schrodinger equation in a diabatic framework. The probability of desorption which is the predicted 
2PACC signal has been calculated as a function of the pulse parameters. 
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I. INTRODUCTION 



Controlling the outcome of a chemical reaction has been the ultimate goal of chemistry. 
Coherent control, a new addition to the quest, is based on exploitation of quantum inter- 
ferences of matter-waves to build a constructive interference in the desired channel and a 
destructive interference in all other channels. To carry out such a task the molecular matter- 
wave has to exhibit a coherent property. To date, all experimental applications of coherent 
control have been based on imprinting the coherent properties of a light/optical source 
onto the matter to be controlled I U [M. With the experimental realization of Bose-Einstein 

nnnn 

condensation |3J, IJ. |5|. |(| . a new source of coherent matter-waves has become potentially avail- 
able. The present paper explores the direct employment of a matter-wave coherent source 
to achieve the goal of coherent control of a chemical reaction. 

The utilization of a matter-wave source depends on the experimental ability to direct and 
shape the matter- waves from a Bose-Einstein condensation (BEC). The waves in a BEC trap 
constitute a single many-body matter-wave function which is an analogue of a photon field 



in an optical cavity. An output coupler termed atom laser transforms t 



re condensate into a 
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. We propose to 



source of either pulsed or continuous coherent matter [7!. S. licl 
employ this matter-wave for a surface mediated chemical reaction carried out by a two-pulse 
atomic coherent control (2PACCj spectroscopy. The present paper elaborates on previous 
preliminary studies (See Refs. |l4|, [if 



BEC has been accomplished for a growinglist of atoms { Cs (Ref. Ilfi), Rb (Ref.lSj), 
85 Rb (Ref. 17), 23 Na (Ref. \fL 7 Li (Refs. flflQ, 6 Li (Ref. Q), K (Ref. 19), H (Ref. Q 
and meta-stable He (Refs. |2l|, |22)}. All of the above condensates and additional systems 
to be prepared, are candidates for a source for a coherently matter-wave. 

The principle of interferin g op tical waves in a molecular framework has neatly been 



demonstrated by Scherer et al. 



23, 



24( . Using a sequence of two optical pulses, population is 



transferred from the ground to the excited electronic state in the iodine molecule. The first 
pulse promotes a portion of the wave function to the excited state, which then evolves in 
time. The second pulse, which is phase-locked relative to the first one, excites an additional 
wave packet to the excited state.. The intramolecular superposition of the two wave packets 
is subject to quantum interferences, either constructive or destructive. This gives rise to 
larger or smaller population on the excited state for a given time-delay between the pulses. 
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The outcome is then controlled by the optica. 



the evolution of the initial wave packet 



23, 



24 



phase difference between the two pulses and 
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26]. 



The objective of the study is to control a bimolecular surface reaction using a coherent 
source of matter-waves shaped as two phase-locked pulses. The target substrate IS cl sticl - 
tionary atom adsorbed on a surface. The outcome of the reaction is a desorbed molecule in 
the gas phase composed from one atom from the matter-wave pulse and one surface atom. 
The method is therefore termed two-pulse atomic coherent control (2PACC) spectroscopy. 
The yield of the reaction is controlled by the time delay between the two atomic pulses and 
the relative phase between them. 

The majority of studies in coherent control have concentrated on unimolecular processes. 
An exception is the study by Abrashkevish et al 27] which theoretically considered a 3D- 
atom-diatom reactive scattering process. The 2PACC spectroscopy is aimed at the more 
difficult task of controlling a bimolecular reaction described by A + B — > C . 

In contrast to a two-pulse photon field the matter waves are dispersive. This means 
that the source of the matter-waves, the BEC condensate, has to be positioned ver y cl ose 



to the target surface. Such a device is realized in the so-called atom chips[2i 



ver y cK 
0,0 



or 



surface micro-traps [31, 32]. In the atom chips for example the BEC has been placed a 



few hundred microns above a metal surface 



291 ] . The additional step required to realize the 



current proposal is a two-pulsed output coupler directing a coherent pulse of atoms toward 
the surface. 

A coherent optical spectroscopy employed in surface science is two-photo emission (2PPE) 
which has served as the inspiration for the present proposed 2PACC. In 2PPE, two photons 
from a laser are applied to a surface with or without adsorbates. In 2PPE, the first optical 
pulse excites a surface electron to an intermediate level of either an image potential state or 
an unoccupied molecular orbital of the adsorbate. The second photon exploits this excitation 
to create a measurable outcome. The method has been used to characterize the energetics 
and dynamics of electrons, atoms and molecules on surfaces 0, Q- S- 3(1 137I. The 
application of two laser pulses allows for probing the nuclear dynamics of the electronic states 
which could not be studied if only one pulse was applied. Recently, Petek et alj^l have 
demonstrated that applying a sequence of phase-locked optical pulses coherently controls 
the motion of an alkali atom on the metal surface and thereby the desorption process. 

The bimolecular surface mediated reaction to be controlled by the 2PACC spectroscopy is 
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the Eley-Rideal reaction. The atomic pulses consist of a coherent atomic source of hydrogen 
or alkali atoms which collide with a hydrogen atom chemisorbed on a Cu(lll)-surface. The 
outcome is an ejected molecule in the gas phase of either H 2 or an alkali hydride. It will be 
shown that the application of a second atomic pulse induces the matter-wave interferences 
that enable enhancement or suppression of the desorbing flux of molecules relative to a single 
atomic beam experiment. The quantum interferences are controlled by the time delay and 
the relative phase between the two atomic pulses. 
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;ical studies of Eley-Rideal reactions have been performed by Jackson et al. 
for an incident hydrogen atom beam interacting with various coated surfaces. 
They have observed that hot atom processes, where the impinging atom becomes trapped 
onto the surface, play an important role. The trapped atoms can react with the adsorbate 
or they can eventually dissipate enough energy through collisions with the adsorbates to 
become immobile. These dissipative forces cause decoherence and therefore will suppress 
the 2PACC control. In this paper the energy and charge transfer between the surface and 
the adsorbates are not included. 

The dynamics of the 2PACC spectroscopy model is explore d by solving the time depen- 
dent Schrodinger equation by a Newtonian propagation method 44|. The matter- wave pulses 
are presented by two phase related Gaussian-like wave-packets. The atomic and molecular 
as well as their interactions with the surface are described in a diabatic framework. 

The methodology of the 2PACC spectroscopy is presented in Section|n] In Section ITTTI the 
dynamics of 2PACC is presented. In Section HVl two simulations of 2PACC spectroscopy are 
studied one with a hydrogen source and another with a lithium atom laser both impinging 
on a Cu-surface with chemisorbed hydrogen atoms. Conclusions and a general outlook are 
presented in Section IVl 

II. THE MODEL 

The target of control is the Eley-Rideal reaction described by 

Y + H/Cu(lll) ->• YH + Cu(lll) . (1) 

The atom Y is from the atom laser source, which in the present model consists of either 
hydrogen or alkali atoms. The matter-wave of Y is directed to a Cu(lll)-surface with low 
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coverage chemisorbed hydrogen atoms. When the wave function of Y overlaps with that 
of the adsorbed atom, interaction is expected, leading to a recombination that forms the 
YH molecule. If the newly formed molecule has sufficient energy it will eventually desorb 
from the surface to the gas phase. By shaping the wave function of Y, the yield of the 
desorbing molecules becomes controlled by constructive or destructive interference. The 
simplest controlled wave function is obtained by a sequence of two pulses where the time 
delay and relative phase are the control variables. It is anticipated that atom lasers will 
be able to produce such coherent matter-wave pulses by an output coupler of a coherent 
source - the Bose-Einstein condensate (BEC). For example a coherent atomic pulse can be 
realized by using an optical laser source to transfer the atom from a trapped (BEC) to an 
untrapped state (the atom laser). The relative phase between the two optical pulses can be 
utilized to phase-lock the two pulses of the atom matter-wave. A time delay between the 
two optical laser determines also the time delay between the two pulsed atomic laser. The 
coherent properties of interest are projected onto this wave function. In this study the wave 
function of an atom laser composed of a single atom is represented by two Gaussian wave 
functions with a well defined time-delay and phase. 

The binding energy for a hydrogen atom on a Cu(lll)-surface is roughly 2.4 eV. Since 
the bond energy of H2 is approximately twice this value, the reaction is very exothermic. 
For this reason the desorbing molecule is expected to be vibrationally hot due to the large 
exothermicity. As one proceeds through the list of atoms (Y={H,Li,Na,K,Rb,Cs}) the 
reaction changes from exothermic to endothermic. If the Cu(lll)-surface is replaced with 
an Al(110)-surface all the reactions becom e ag ain exothermic since the binding energy of 
the hydrogen on this surface is only 1.8 eV |45j ]. 

The present modeling of the Eley-Rideal reaction restricts the motion of the two atoms 
to a collinear configuration normal to the surface. We consider the reaction between an 
impinging atom of mass m y located at z y and a target atom of mass located at Zh 
which is initially adsorbed on a flat static surface. The coordinate set (zh, z y ) is transformed 
to a new set of coordinates (r, Z) where r = z y — Zh is the intramolecular distance and 
Z = (rrihZh + m y z y )/M is the center of mass coordinate. The total and reduced mass of 
the system are denoted by M and fi, respectively. Effects arising from phonons and surface 
corrugations are not included in this model. 
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A. Two-dimensional potential energy surfaces 



To understand the dynamics of either a dissociation reaction of a diatomic or a recombina- 
tion reaction of two atomic species on a surface, potential energy surfaces (PES) representing 
the physics of breaking and forming a chemical bond is required. A comprehensive overview 
of the historical development of potential energy surfaces is given in Ref. Q- 

An electronic structure model for the dissociation/recombination of the H 2 /Cu system 



47, l48l. Their model is based on a small 



has been suggested by Holloway and co-workers 
complex of the form Cu 2 H 2 in a planar C% v geometry. The intramolecular coordinates are 
r which is the distance between the two hydrogen atoms, and Z the distance between the 
midpoint of H-H and Cu-Cu. For large value of Z there is no interaction between the 
two units, H-H and Cu-Cu, meaning that each unit has a separate set of fully occupied 
molecular orbitals. As Z decreases the two units approach each other, the orbitals with 
similar symmetry with respect to the bond center mix, e.g. the parallel-bonding orbital of 
H-H and Cu-Cu interfere with the parallel anti-bonding of Cu-H. Due to these interactions 
the orbital energies shift. This causes a crossing between the "parallel" and "perpendicular" 
orbitals with the result that the H-Cu bonding becomes energetically favorable. 

The PES of the reactant surface has been constructed from a Morse potential describing 
the hydrogen- hydrogen bond and a repulsive potential as the hydrogen molecule approaches 
the surface. In the product channel the two separated hydrogen atoms are chemisorbed on 
the surface, this bond is described by a Morse potential. The hydrogen-hydrogen repulsion 
on the surface is represented by an exponential repulsive potential. 



In the diabatic representation, the potential is represented by a 2x2 V-matrix 



V(r,Z) 



V RR (r,Z) V RP (r,Z) 
V PR (r,Z) V PP (r,Z) 



(2) 



The diagonal elements are the potential energy surfaces of the reactant and the product 
states and the off-diagonal elements is the non-adiabatic couplings between them. Upon 
diagonalization of the 2x2 V-matrix two new adiabatic PES are obtained. In the previous 
studies of Eley-Rideal reactions the adiabatic PES has been used. We will briefly discuss 
the idea behind them. 

4^1 used the London-Eyring-Polanyi-Sato (LEPS) potential energy 



Jackson et al. 
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surface for modeling the Eley-Rideal reactions. The LEPS potential is given by 
V = U m + U a + U b ± y/Al+[A a + A b }2-A m [A a -A b ] , 



(3) 



where U and A are the Coulomb and exchange contributions. The intramolecular contribu- 
tions are denoted U m and A m , whereas the other terms describe the atomic interactions with 
the surface. Here, only the lowest adiabatic PES are considered. The PES of the molecule 
as well as the interactions between the surface and a single atom have been calculated by 
Density Functional Theory (DFT). They have been fitted to the functional forms, U and 
A, which are based on Morse-like attractive and repulsive terms which decay exponentially 
with increasing distances. Surface effects have been introduced through expanding the Morse 
parameters representing the surface in a Fourier series of the reciprocal lattice vectors. 

The LEPS-potential is a member of a class of diatomic in molecules (DIM) potential 
forms. For surface reactions, an extension called embedded diatomic in molecules (EDIM) 

and recently 

used for the description of the motion of N 2 on a ruthenium-surface {5 1| . In the EDIM the 
intramolecular interactions are modeled by Morse and anti- Morse potentials which represent 
the singlet and triplet electronic state of the diatomic molecule in the gas phase. The 
interaction between the atoms in the gas phase and the surface is modeled by embedded 
atom model (EAM), in which the atom experiences an average charge density from the 
surface. 

In the 2PACC spectroscopy the coherent properties of the encounter is intricate therefore 
the Born-Oppenheimer approximation in which the reaction takes place on a single potential 
energy surface has to be replaced with a non-adiabatic framework. In the 2PACC dynamics 
a diabatic frame is chosen consisting of two potential energy surfaces. Specifically, for the 
Eley-Rideal reaction the diabatic PES are constructed from two atomic or one molecular 
chemical species interacting with the surface. In the reactant channel the interaction between 
the stationary hydrogen atom and the surface is strong due to a chemical bond. We are using 
a semi-empirical functional form for the PES. The impinging atom experiences a repulsive 
force from the adsorbed atom as well as a long-range attraction from the surface due to 
polarization forces. In the product channel the newly formed molecule is attached to the 
surface by a weak bond induced by the polarization forces between the surface and the 
molecule. The PES for the reactant and product channels have the following functional 
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forms 

V RR (r, Z) = V£ M + V* M + V? H (4) 
V PP (r, Z) = VVVm + Vf H , (5) 

where the index M represents the interactions with the surface. The last terms of Eqs.flHEJ) 
represent the intramolecular interaction between the two gas phase atoms, H and Y. The 
other terms represent the interaction between an atom or a molecule with the surface. In 
the following section we will discuss each of the individual terms in the reactant and product 
PES. 

TABLE I: The parameters for PES for the Eley-Rideal reaction on a Cu(lll)-surface. The 
hydrogen atom is chemisorbed on a hollow site of the surface. 

hydrogen molecule £^#=4.505 eV a Z e =2.0A A = 0.2 

rf=0.74lA a rf=0.754A 6 
a 9 HH =2.2k- 1 a^ H =2.1lA- 1 

chemisorption D HM =2.2>M(N a ^ M =0.916A a a HM =l-7^~ l a 

Physisorption A M /eV &m/A _1 C M /eVA~ 3 

H 2 or H 600 3.8 3.5 

LiH or Li 650 3.15 7.0 

NaH or Na 760 2.7 12.5 

KH or K 850 2.4 20.0 

RbH or Rb 930 2.23 27.0 

CsH or Cs 950 2.1 34.0 
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1. Intramolecular interactions 



Asymptotically, the intramolecular interaction potentials, Vyn an d Vyh m the reactant 
and product channels become the triplet and singlet electronic states of the molecule, YH, 
when both atoms are far from the metal surface. The PES of the alkali hydride molecule 
was based on the intramolecular gas-phase singlet and triplet electronic states which have 
been evaluated by a multi-configuration self-consistent field calculation by Geum et al.[53]. 
These potentials have been fitted to a cubic spline interpolation function. 

For the hydrogen molecule the singlet and triplet PES are given by a Morse and an 
anti-Morse potential 

Vg H (r,Z) = D HH ([1 - e {-<*HH(Z)(r-r« HH (Z))}}2 _ ^ 

Vg H (r,Z) = ^^([l + ^^f-l) . 

The dissociation energy of the hydrogen molecule is denoted D HH and (1 — A)(l + A) is 
the Sato factor. The equilibrium distance between the two hydrogen atoms, r e HH , and the 
coupling strength, (Xhh, have been obtained as a function of Z by an interpolation between 
values for the gas phase to ones for the adsorbate: 

a HH (Z) = a a iu for Z < Z e 

a HH {Z) = a g HH - {a 9 HH - a a x H )T 4 (Z, Z e , 2b HM ) for Z > Z e (6) 

and 

r% H (Z) = r* - (rf - rf)T,(Z, Z e , 2b HM ) . (7) 

The superscripts g and ad indicate the gas phase and the adsorbed state. Table U gives the 
parameters used for the PES for a Cu(lll)-surface covered by a hydrogen atom adsorbed 
on the hollow site. 

2. Surface interactions 

In the product channel a molecule (YH) is formed and eventually transferred to the gas 
phase. The molecule is coupled to the surface by a weak bond induced by the polarization 
forces. The interaction potential, Vy H _ M , is represented as a sum of two terms - a short 
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FIG. 1: The potential energy surface representing the physisorption of a hydrogen or an alkali 
hydride molecule on a Cu(lll)-surface are shown as a function of the distance between the surface 
and the mass center of the diatomic. 

range repulsive term and a weak long-range attraction between the molecule and the surface 

Vyh-m = A M e^ bMZ) " - r 4 (Z, 0, 2b M )) . (8) 

The incomplete Gamma function which turns off the attraction as the Y-atom approaches 
the surface is given by: 

k=m / n£. 

r m (x, so, a) = J2 1l ■ ( Q ) 

k=0 

The parameterization of the physisorption of the molecular interaction with the surface 
has been chosen such that the minimum energy e.g. the physisorption energy is observed 
at Z e ={r e HM + myTyi^lM. The location of the minimum energy is moved further away 
from the surface as the mass of the impinging atom increases. As Y becomes heavier, the 
polarization force between the surface and the alkali atom increases, thereby the binding 
energy of the physisorption bond becomes stronger. The amplitudes of the repulsion, Am 
and attraction, Cm increase as one proceeds down the list of atoms, {H,Li,Na,K,Rb,Cs}. 
The coupling strength, 6m, is inversely proportional to Z e . The parameters are constructed 
so that the physisorption energy of YH increases slightly as one proceeds down the list of 
atoms {H,Li,Na,K,Rb,Cs}. The physisorption potentials are shown in Figure ^ and the 
parameters are given in Table |l] 
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To model the interaction between the impinging atom in the reactant channel a physisorp- 
tion potential has been used. The same set of parameters for the atomic interaction as for 
the molecular physisorption interactions have been chosen but the coordinate, Z, in Eq.fjSJ) 
is replaced by the location of the impinging atom, z y . 

In the reactant channel the interaction between the chemisorbed atom and the surface is 
described by a Morse potential 

Vg M = D HM {[1 - e (-HM(^-4«))]2 _ i) ; ( io) 

where D HM is the dissociation energy of the atom-metal bond and z e HM is the equilibrium 
bond length. The parameters are given in Table H] for a hydrogen atom adsorbed at the 
hollow site of a Cu(lll)-surface. 

B. Wave function of the atomic pulse 

Initially, when the impinging atom is far away from the surface and the adsorbed atom, 
the wave function representing the outcome of the total system is written as a product of a 
wave function on each atom, 

V atom (r, Z) = Nip h (z h , z y = z s ) ^ y {z y ) ; ^ y {z y ) = e -(*»-*)V*v (11) 

where N is a normalization factor. The wave function of the adsorbed atom, iph, was chosen 
as the lowest energy eigenstate of the PES of the reactant channel along the one-dimensional 
path with a fixed distance, z s , between the approaching atom and the surface. Considering 
only the lowest vibrational states would be sufficient since in a given experimental situation 
one would expect that the temperature of the surface is ultracold. The wave function of 
the approaching atom, ip y , is represented by a Gaussian wave function centered at z y = z s 
and with a variance, a y . The momentum of the approaching atom is denoted by k y = 
±y^2m y E/h 2 which is related to the energy E of the propagating atom laser. The sign of the 
momentum operator determines the orientation of the atom laser: For negative momentum, 
the atom moves toward the surface. The variance is related to the dispersion of the atom 
laser. 
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III. 2PACC DYNAMICS 



The 2PACC spectroscopy of an Eley-Rideal reaction consists of the dynamics on the 
coupled reactant and product potential energy surfaces. The wave function of the system is 
described by the vector 

'* R (r,Z,t) 
* P (r, Z,t) 

where the wave functions in the reactant and product channels are denoted by the index 
{R, P}. The Hamiltonian of the system is represented by 



*(r, Z,t) 



(12) 



+ 



H(r, Z) = T{r, Z) + V(r, Z) 

T RR (r,Z) 
T PP (r,Z) 

The diagonal elements of the Hamiltonian have the following form 

h 2 d 2 h 2 d 2 



V RR (r,Z) V RP (r,Z) 
V PR (r,Z) V PP (r,Z) 



(13) 



H u (r, Z) = T«(r, Z) + V u (r, Z) = -- 



V«{r,Z) 



{R,P} (14) 



2M dZ 2 2/x dr 2 

where the first two terms are the kinetic energy operator for the two degrees of freedoms (r 
and Z) and the last term is the potential energy function which was described in Section 
III Al The off-diagonal elements of the Hamiltonian represent the non-adiabatic coupling 
between the two channels which are described by 



V RP (r,Z) 



V PR (r,Z) 



pz 



(15) 



where (3 is the non-adiabatic coupling strength and ry H is the equilibrium distance between 
the two atoms in the singlet molecular PES. This representation ensures that the electron 
density of the metal decays exponentially outside the surface into the vacuum. It is important 
to note that the adiabatic PES of the Eley-Rideal reaction which have been used by Jackson 



et al. 



3, 41, 42, 3 



can be obtained by diagonalization of this 2x2 V-matrix. 
The dynamics of the 2PACC was followed by solving the time- dependent two-channel 
Schrodinger equation which is given by 

9* 



%h- 



dt 



(16) 



In the two-pulse atomic coherent control (2PACC) spectroscopy the first atom laser rep- 
resented by the wave packet (Eq.([TTj)) is initialized at the reactant PES at a time t = e.g 
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^atom(r,Z) \ ^ 

/ 

This initial wave packet evolves in time and after a time delay, At, the second atom pulse 
is introduced given by 

^atom,2(r, Z) = ^ a tom(r, Z) exp (-Z0) (18) 

where 9 £ [— 7r; 7r] is the relative phase between the two atomic pulses that describes the 
coherence between them. As the wave packet propagates population is transferred from the 
reactant to the product PES through the non-adiabatic coupling. 



*(r,Z,0) 



A. Computational Method 

The wave function is represented on a two-dimensional grid. First, the wave function 
of the chemisorbed hydrogen was calculated. Using an imaginary time propagation [54] the 
one-dimensional wave function of the hydrogen atom, ipjj along the line with a fixed distance 
between the approaching atom and the surface, z y =6.8A has been relaxed to its vibrational 
ground state. Representing the wave function of the impinging atom as a Gaussian wave 
function, the total wave packet Eq. ljll)) is then initialized. 

The dynamics of the 2PACC is obtained by propagating the initial wave function by 
exp (— iHt)ip(0), in which the time-evolution operator exp(-iHt) is expanded by New- 
tonian interpolation polynomials with Chebychev sampling points 55j . The kinetic energy 
operator has been evaluated using the fast Fourier transformation technique jo^i Isfij ]. The 
parameters used in the wave packet propagation are displayed in Table |HI 

The scattered wave function is removed at large values of r and Z by complex absorbing 



potentials 



571 ] which prevent reflection and transmission at the end of the grid. The overall 



potential can be written as 

V(r,Z) = V(r,Z) + V abs (r) + V abs (Z) if r-A r <r<r max and Z - A z < Z < Z max 
= V(r,Z) + V abs (r) if r - A r < r < r max and Z < Z - A z 

= V(r, Z) + V abs (Z) if r <r - A r and Z - A z < Z < Z max 

= V(r,Z) if r<r-A r and Z < Z - A z (19) 
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TABLE II: Computational parameters for the wave packet propagation of the 2PACC dynamics. 



time step 
propagation time 
time steps 



5t=0.097 fs 



tmax — 484 fs 



N t =5000 



grid points along r 
grid spacing along r 
grid starts at 
grid points along Z 
grid spacing along Z 
grid starts at 



N r =256 
Ar=0.0529 A 
rw=0.0529 A 
N z =256 
AZ=0.0529 A 
Z mir) =0.0529 A 



variance 
initial position 
momentum 



<7„=0.280 A 2 
z s =6.82 A 
fc=9.45 A" 1 



absorbing potential 



A r =A z =1.32 A 
y =0.00027 eV 



non-adiabatic coupling 
dividing flux line 



/?=0.027 eV 



■'flux' 



=5.24 A 



where the complex potential is given by V a b s (r) = iv Q (r — (r max — A r )) 2 . The same functional 
form is used for Z. r max is the last grid point and A r is the interval where the complex 
potential is applied. 

As the wave function is evolving from the entry/reactant channel, the non-adiabatic 
coupling term is responsible for transferring the amplitude between the two diabatic surfaces. 
Eventually the molecule desorbs from the metal surface. This leads to an outgoing flux in the 
exit/product channel. The probability of a YH- molecule to escape from the metal surface is 
obtained via a flux-resolved analysis carried out at an asymptotic value of Z=Z 'fi ux —5. 24A. 
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The total accumulated desorbing flux which is the predicted 2PACC signal is computed by 



Nt 



i 

5t Nt 
= — ) Im 



(20) 



dr &*p(r,Zf lux ,ti 



8Z 



where the derivative is evaluated with a Fourier Transform. The integrated flux has been 
determined as a function of the time delay, At, and the phase-relation, 9, between the two 
atom laser pulses. 

Furthermore, the accumulated flux current can be evaluated for each of the vibrational 



states along the dividing line 



Ni 



(21) 



6t 



N t 



£)lm ** n (Z flux ,t) 



d* n (Z,t) 



z=z 



flu 



M ^" 'V " ' dZ 

i 

where j n is the probability current for the wave packet to go into the n'th vibrational state. 
Here, *Sf n (Z,t) is the projection of the wave function onto the vibrational eigenstates, Xn, 
along the dividing line 

* n (Z,t) = J drx* n (r)* P (r,Z,t) . (22) 

The vibrational eigenstatates have been calculated by imaginary time propagation. (See 
Appendix 1X1 for further details). 



IV. RESULTS AND DISCUSSION 



The 2PACC of a hydrogen source is compared to that of a lithium atom laser source im- 
pinging on a surface with chemisorbed hydrogen atoms. The evolution of the wave packet on 
the reactant and the product channel are presented. The 2PACC signal and the vibrational 
analysis of the desorbing molecule are displayed. The control parameters are the time delay 
and relative phase between the two atom laser pulses. 



A. H+H/Cu(lll) 



The dynamics of the 2PACC spectroscopy are demonstrated for a hydrogen atom laser 
source impinging on a Cu(lll)-surface coated with hydrogen atoms. Compared to earlier 
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FIG. 2: The evolution of the wave packet (|vl/#| 2 ) from an atom laser source applied to a Cu(lll)- 
surface coated with hydrogen atoms superimposed on the reactant potential energy surface. Snap- 
shots are shown for a time delay of At=4.84 fs and the phase, 6=— ^rr. The contour values of 
potential energy surfaces are -3, -2.5, -2, -1.5, -1 and -0.5 eV. The PES of the reactant and the 
product channels are shown with blue and black lines, respectively. 
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FIG. 3: The evolution of the wave packet on the product channel (|\l/p| 2 ) are shown for time 
delay of Ai=4.84 fs and phase, 6=— ^ir. Notice the exciting wave packet on the reactant channel 
at 145.1 fs. Parameters used for this calculation are given in Figure [21 



calculations 
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the grid spacing in both degrees of freedom (r,Z) was reduced by a factor 



of 2. The initial starting position of the wave function was z s =Q.8A compared to 13. 6A in 
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Phase [9] 

FIG. 4: The 2PACC signal of a hydrogen atom laser impinging on a surface with chemisorbed 
hydrogen atoms relative to the outcome for one pulse (in %). The signal is shown as a function 
of the time delay (in fs) and relative phase (9) between the two atomic pulses for an initial wave 
packet with energy -2.3 eV. The predicted 2PACC signal is calculated from the integrated flux 
along a dividing line on the product channel at t = 484 fs. Areas of enhancement are in red and 
those of depletion in blue. 

the earlier calculations. At the shorter distance the interactions between the impinging atom 
and the adsorbate and the surface is still negligible. 

In Figures 121 and |3] snapshots of the evolution of the wave packet on the reactant and 
the product surfaces are shown for 6=—^tt and A£=4.84fs. First, the initial wave packet 
(Eq.([lip) is generated and this atom laser pulse evolves in time. The energy of the initial 
wave packet is -2.3 eV and the impinging atom has a kinetic energy of 0.2 eV. After the 
specified time delay a wave packet - the second atom laser pulse - is placed in the original 
position of the first wave-packet (see second snapshot of Figure EJ. The relative phase 
between the two wave packets at the initial position may be different e.g. 8^0. The wave 
function which resembles the two atom laser pulses propagates toward the non-adiabatic 
region where the reactant and the product surfaces intersect. As the wave packet enters 
the non-adiabatic region a part of the wave function is transferred to the molecular state 
where it can exit. This means that the approaching hydrogen atom has reacted with the 
adsorbed hydrogen and a H2 molecule is formed. The non-reactive part of the wave function 
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FIG. 5: The wave packet after the second atom laser pulse has been applied at t=At=14.5 fs for 
different relative phases between the two atomic pulses. 

collides with the adsorbed hydrogen atoms. Either it scatters back to the gas phase or it 
gets trapped on the surface due to the weak polarization forces between the impinging atom 
and the surface (physisorption). 

The population which is transferred through the non-adiabatic coupling to the product 
channel builds up slowly as a function of time. In the beginning the newly formed hydrogen 
molecule is physisorbed on the surface. That is, it is trapped in the potential well. Later 
a fraction of the wave function exits the channel and the hydrogen molecule desorbs from 
the surface to be detected in the gas phase. The desorbing molecule is vibrationally excited. 
The node structure of the exiting wave function shows that the hydrogen molecule desorbs in 
the fifth vibrational state (V=4). This observation will be verified by a vibrational analysis. 

The 2PACC signal, which is the difference between a one pulse and a two pulse desorbed 
molecular outcome is calculated from the total integrated flux at t = 484 fs. Figure |U shows 
the 2PACC signal as a function of time delay and relative phase between the two pulses. 
The amount of control of the 2PACC signal demonstrated in this model is an enhancement 
of 350 % for 6=-\n and At=4.84 fs relative to one atomic pulse compared with a suppression 
of 95% at 9=n and At=9.7 fs. 

The 2PACC signal shows a variation with respect to both control parameters, meaning 
that the outcome of the Eley-Rideal reaction is coherently controlled by the time delay and 
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FIG. 6: The accumulated flux in each of the vibrational states for the one pulse experiment are 
shown for hydrogen (upper figure) and lithium (lower figure). 

the relative phase between the two atomic pulses. The application of the second pulse creates 
a quantum interference with the first pulse. Such interferences are either constructive or 
destructive resulting in increasing or decreasing the flux of the desorbing molecules in the 
product channel. The effect of these interferences can be visualized by considering the wave 
function just after the second atom pulse has been applied. In Figure|S]the total wave packet 
is shown for i=Ai=14.5fs for different relative phases between the two atom laser pulses. 
An elongated structure of the wave packet immediately after the second pulse gives rise to 
destructive interference and decrease of the desorption yield whereas a "node-like" structure 
gives rise to constructive interference and an enhancement of the yield. 

In FigurelHlthe accumulated flux in each vibrational state of desorbing hydrogen molecule 
is shown for a one pulse experiment. The hydrogen molecule leaves the surface predominantly 
in the fifth lowest vibrational excited state (z/=4). This observation agrees with the node 
structure of the leaving wave function on the product channel (See Figure EJ). Jackson 
et alj3] also observed a vibrational distribution of desorbing hydrogen molecules centered 
around z/=4 for an Eley-Rideal reaction of two hydrogen atoms on a Cu(lll)-surface in a 
collinear quantum mechanical calculation. 

The extent of control in the accumulated vibrational flux between the two and one pulse 
scenario is shown in Figure For all the calculations the vibrational distribution is centered 
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FIG. 7: The percentage deviation in the accumulated flux in the vibrational states from the one 
pulse experiment are shown for four different relative phases (9) between the two pulses. For each 
of the phases the accumulated fluxes are displayed for four time delays 1) Ai=4.84 fs (dark blue), 
2) Ai=9.68 fs (light blue), 3) At=14.52 fs (yellow) and 4) At= 19.35 fs (red). The accumulated 
flux in each vibrational state for a one pulse experiment is shown in Figure H3 

around z/=4. An enhancement of flux in the lowest vibrational states of the hydrogen 
molecule (z^<3) is correlated with an increases of the total 2PACC yield. On the other hand 
an enhancement of the higher lying vibrational states (z^>4) is anticorrelated with the total 
2PACC yield. 



20 



12345678 12345678 




12345678 12345678 
r in Angstrom r in Angstrom 



FIG. 8: Evolution of the wave packet on the reactant surface for the lithium atom laser impinging 
on a H/Cu-surface. The time delay and the relative phase between the two atom laser pulses are 
A£=9.7fs and 9=-^ir. The contours given for -3, -2.5, -2, -1.5, -1 and -0.5 eV for the potential 
energy surfaces of the product and the reactant channels are shown with blue and black lines, 
respectively. 

B. Li+H/Cu(lll) 

Next, the lithium with an adsorb hydrogen Eley-Rideal reaction on a Cu(lll)-surface 
is studied. A smaller probability of desorption for lithium hydride molecule is expected. 
This is due to the reduced energy difference between the product and the reactant and a 
stronger physisorption interaction between the surface and the lithium hydride. For a single 
matter-wave pulse, the desorbing flux of lithium hydride is suppressed by 10 -9 compared 
to desorption of hydrogen. The total energy of the initial wave packet is -2.4 eV where the 
initial kinetic energy of the impinging atom is 0.027 eV. 

The evolution of the wave packet on the reactant and the product PES is shown in Figures 
|H]and|ni The dynamics of the 2PACC spectroscopy with a lithium atom laser is slower than 
the one for a hydrogen atom laser since the mass of the impinging atom increases by almost 
a factor of 7. As the wave packet on the reactant channel approaches the metal surface, 
part of the population is transferred to the product channel. The magnitude of population 
transfer is much smaller for lithium than for hydrogen as expected. The part of wave packet 
on the reactant channel which does not react is primarily trapped on the surface due to 
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FIG. 9: Evolution of the wave packet on the product PES for lithium atom laser impinging on 
a H/Cu-surface. The time delay and the relative phase between the two atom laser pulses are 
At=9.7fs and Q=-\k. 

polarization forces and only a very small part of the wave function scatters back to the gas 
phase. The opposite was observed for the hydrogen case. 

The predicted 2PACC signal is shown in Figure ITU1 as a function of time delay and phase 
between the two pulses. For certain values of time delay and phase between the two pulses 
the second pulse induces a large constructive interferences. These interferences enhance the 
probability of desorption for LiH by 2100% relative to a single-pulse experiment. 

The desorbing molecule, LiH, is vibrationally excited. This is because the bond length 
LiH is considerably longer than chemisorption bond between the hydrogen and the Cu(lll)- 
surface. Along the dividing line where the desorbing flux is collected, the fourteen lowest 
vibrational eigenstates of the lithium hydride have been evaluated. For a single pulse ex- 
periment (See Figure a broad distribution of the vibrational states with a maximum at 
u=8 is observed. In Figure HP the enhancement /suppression in the integrated vibrational 
flux relative to the single-pulse scenario is shown. The deviation in the integrated flux as 
a function of the phase reflects the 2PACC signal. If there is a large enhancement of the 
2PACC signal the accumulated flux in each vibrational state increases. 
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FIG. 10: The 2PACC signal. The enhancement relative to the output of a single pulse for lithium 
atom laser applied to a H/Cu(lll)-surface. The control parameters are the time delay and phase 
between the two pulse. The predicted 2PACC signal is calculated from the integrated flux along a 
divided line on the product channel at t = 484 fs. Areas of enhancement are in red and depletion 
are in blue. 

V. SUMMARY AND CONCLUSIONS 

A two-pulse atomic coherent control (2PACC) spectroscopy has been presented in this 
study. The coherent properties of an atom laser have been used to control a surface mediated 
chemical reaction, A + B — > C , where the A wave function is generated from an atom laser 
and the target atom B is immobilized by the surface. Two atom laser pulses are applied 
creating an interference pattern with each other. These interferences are the essence of the 
coherent control of the chemical reaction. The two pulses are necessary since the initial 
relative phase between the atom A and atom B on the surface is arbitrary j^. 0]. The 
control knobs used were the time delay and a relative phase between two atom laser pulses. 
A significant enhancement of the 2PACC signal relative to single atom laser pulse was 
obtained. 

The Eley-Rideal reaction forming H 2 is highly exothermic whereas the recombination 
reaction leading to LiH is almost thermoneutral. Using a coherent matter-wave source 
enables us to coherently control the recombination reaction. The 2PACC spectroscopy 
opens up reaction channels which are almost closed to incoherent sources. The reaction 
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FIG. 11: The percentage deviation in the accumulated flux in the vibrational states relative to the 
single pulse simulation. Four different relative phases (6) between the two pulses are shown. For 
each of the phases the accumulated fluxes are displayed for four time delays 1) Ai=38.7 fs (dark 
blue), 2) At=43.5 fs (light blue), 3) At=48.5 fs (yellow) and 4) Ai=53.2 fs (red). The accumulated 
flux in each vibrational state for the single pulse simulation is shown in Figure H3 

yield by Eley-Rideal reaction of lithium atom laser with adsorbed hydrogen atom has been 
enhanced by more than 2100% compared to a one-pulse reference. Here, we manipulated the 
wave function of one atom transferred from a trapped state in the BEC to an untrapped and 
propagating state. The coherent properties of interest are projected onto this wave function. 
In this study a double Gaussians wave function was used where the source of coherent control 
is the time delay and relative phase between the two Gaussian wave functions. The design of 
the output coupler is the experimental "bottleneck" in coherent control with matter- waves. 
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In analogy with the optical pulse shaper (fulfill, a more elaborated matter-wave pulse shaper 
controlling the amplitude and phase of the matter-wave would eventually lead to a reaction 
probability of one. 

The design of the wave function could further be used to direct the leaving molecule 
into a specific rovibrational state. The 2PACC spectroscopy could also be employed for 
endothermic reaction such as Eley-Rideal reaction between a adsorbed hydrogen atom and 
sodium (rubidium, cesium). 

Additional coherent manipulations could be possible if the atomic matter-wave pulse is 
replaced by a molecular pulse where the source is a molecular Bose-Einstein condensate. 
For a heterogeneous molecular matter-wave the branching ratio between the two reaction 
channels could be controlled. From the prospective of coherent control methodology the 
current study demonstrates control of a recombination reaction. This is in contrast to optical 
coherent control which has been almost exclusively limited to the control of unimolecular 
reactions such as photo-dissociation. 

Control of a recombination reaction is another direction where both atoms come from the 
same BEC source. In this case their initial relative phase is well defined while the surface 
serves to break the symmetry. Since both sources are the coherent the reaction product is 
also coherent. The ultimate goal is to design the wave functions of the atom lasers such 
that the reaction probability is unity for one specific rovibrational states of the product. 
This would lead to the formation of a stable molecular BEC. Control on these lines is under 
investigation. 
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APPENDIX A: VIBRATIONAL EIGENSTATES 



Vibrational eigenstates, Xn, are solutions of the one-dimensional, one channel time- 
independent Schrodinger equation for a fixed distance, Zf [ux , from the surface 

H PP (r; Z f i ux )xn(r; Z flux ) = E n (Z flux )xn(r; Z flux ) . (Al) 

Asymptotically when Zfi ux — > oo, the eigenstates represent the vibrational states of the 
diatomic molecule in the gas phase. The method for evaluating the eigenfunction is based 
on propagating a trial wave function accor ding to the time dependent Schrodinger equation 
in imaginary time t — it. A Gaussian filter |§2| has been used to project the eigenstates out 
in an energy region around an initial guess, e, 

Xtnai(r, r) = e-^ H -^ hAE Xtnal (r, 0) . (A2) 

The parameter r has the dimension of time while AE is the energy range covered by the 
numerical procedure which appears as a scaling factor for the normalization of the Hamilto- 
nian in the Newtonian time propagation. This procedure can be interpreted as a relaxation 
of the ground stated of a modified Hamiltonian H* = 4(H — e) 2 /AE. Convergence onto a 
specific target eigenstate can be estimated from the energy dispersion relation 

D(r) = y/( X (r) | if 2 | X (r)> - <x(r) | H | X (r)) 2 (A3) 

which decreases uniformly by increasing the "purity" of an eigenfunction and vanishes in 
the limit of an exact eigenstate. If the energy dispersion, D{r) is less than 1.0E-6 the 
propagation is stopped and the wave function is a vibrational eigenstate of the diatomic 
molecule, Xn for a fixed distance to the surface. 
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